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Besides their astrophysical interest, compact stars also provide an arena for understanding the 
properties of theories of gravity that differ from Einstein’s general relativity. Numerous studies have 
shown that different modified theories of gravity can modify the bulk properties (such as mass and 
radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed 
though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. 
Motivated by this observation, in this paper we take an alternative route and construct a stellar 
structure formalism which, without adhering to any particular theory of gravity, describes in a simple 
parametrized form the departure from compact stars in general relativity. This “post-Tolman- 
Oppenheimer-Volkoff (TOV)” formalism for spherical static stars is inspired by the well-known 
parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable 
correction terms to the fully relativistic TOV equations. We show how neutron star properties are 
modified within our formalism, paying special attention to the effect of each correction term. We 
also show that the formalism is equivalent to general relativity with an “effective” (gravity-modified) 
equation of state. 

PACS numbers: 04.40.Dg, 04.50.Kd, 04.80.Cc, 04.25.Nx, 97.60.Jd 


I. INTRODUCTION 

Neutron stars play a special role among astrophysical 
objects, because they are excellent laboratories for mat¬ 
ter under extreme conditions (unlike black holes) and 
also excellent laboratories to probe strong gravity (un¬ 
like ordinary stars or white dwarfs) [1] . For these reasons 
neutron stars are among the main targets of future ob¬ 
servatories, such as SKA [2], NIGER [3], LOFT [4] and 
AXTAR [5]. These experiments have the potential to 
measure neutron star masses and radii to unprecedented 
levels [6-8]. If general relativity (GR) is assumed to be 
the correct theory of gravity, the observed mass-radius 
relation will constrain the equation of state (EOS) of 
matter at supranuclear densities, which is inaccessible 
to laboratory experiments [9-14]. A procedure to re¬ 
construct the EOS from observations of the mass-radius 
relation (working within GR) was developed in a series 
of papers by Lindblom and collaborators [15-17]; see [18] 
for a review. 

Besides their interest for nuclear physics, neutron 
stars are also unique probes of strong-held gravitational 
physics. For any given EOS, theories that modify the 
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strong-held dynamics of GR generally predict bulk ob¬ 
servable properties (neutron star mass, radius, moment 
of inertia and higher multipole moments) that are differ¬ 
ent from those in Einstein’s theory. However, a survey 
of the literature on neutron stars in modihed theories of 
gravity (see e.g. Table 3 of [1]) reveals a high degree of 
degeneracy in the salient properties of relativistic stars. 
As we show in Fig. 1, if we assume a nuclear-physics mo¬ 
tivated EOS (specihcally, EOS APR [19] in the hgure), 
modihcations in the gravity sector are usually equiva¬ 
lent to systematic shifts of the GR mass-radius curves 
towards either higher masses and larger radii (as in the 
case of scalar-tensor theories [22, 24]), lower masses and 
smaller radii (as in the case of Einstein-dilaton-Gauss- 
Bonnet [20, 25] and Lorentz-violating theories [21, 26]) or 
both, as in Eddington-inspired-Born-Infeld gravity with 
different signs of the coupling parameter [23, 27]. 

Systematic shifts in the mass-radius relation could be 
attributed either to the poorly known physics controlling 
the high-density EOS, or to modifications in the theory 
of gravity itself. This EOS/gravity degeneracy is intrin¬ 
sic in all attempts to constrain strong gravity through 
astrophysical observations of neutron stars: chapter 4 
of [1] reviews various proposals to solve this problem, 
e.g. through the recently discovered universal relations 
between the bulk properties of neutron stars [28-31]. 

In any case, different gravitational theories span (at 
least qualitatively) the same parameter space in terms of 
their predictions for relativistic stellar models. Gravity- 
induced modifications usually look like smooth deforma- 
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dormant for more than thirty years. In the intervening 
time the PPN parameters have been extremely well con¬ 
strained by Solar System and binary pulsar observations 
at IPN order (see [49] for a review of current bounds). 

In this paper we build a phenomenological post-TOV 
framework by considering IPN and 2PN order correc¬ 
tions to the TOV equations. Our strategy is, at heart, 
quite simple: from a suitable set of PPN hydrostatic equi¬ 
librium equations we isolate the purely non-GR pieces. 
These PPN terms are subsequently added “by hand” to 
the full general relativistic TOV equations, hence produc¬ 
ing a set of parametrized post-TOV equations (cf. [50] 
for a similar “post-Einsteinian” parametrization in the 
context of gravitational radiation from binary systems). 
The formalism introduces a new set of 2PN parameters 
that are presently unconstrained by weak-field experi¬ 
ments, and that encompass the dominant corrections to 
the bulk properties of neutron stars in GR in a wide class 
of modified gravity theories. 


FIG. 1. The gravity-theory degeneracy problem. The mass- 
radius relations in different modified theories of gravity for 
EOS APR [19]. Masses are measured in solar masses, and 
radii in kilometers. The theory parameters used for this 
plot are: a = 20 Mq and fP = 1 (Einstein-dilaton-Gauss- 
Bonnet [20]); ci 4 = 0.3 (Einstein-Aether [21]); j3 = —4.5 
(scalar-tensor theory [22]) and k = ±0.005 (Eddington- 
inspired-Born-Infeld gravity [23]). Even if the high-density 
EOS were known, it would be hard to distinguish the effects 
of competing theories of gravity on the bulk properties of neu¬ 
tron stars. 


A. Executive summary 


Since this paper is rather technical, we summarize our 
main conclusions here. The core of our proposal is to 
use the following set of “post-TOV” equations of struc¬ 
ture for spherically symmetric stars (from now on we use 
geometrical units G = c = 1): 


dp 

dr 



pm 

^2 


{V1+V2), 


(la) 


tions of the general relativistic predictions. A notable 
exception are cases where nonperturbative effects induce 
phase transitions, as in the “spontaneous scalarization” 
scenario first proposed in [32], where modifications only 
occur in a specific range for the central density. 

With the possible exception of nonperturbative phase 
transitions, these considerations suggest that the broad 
features of a large class of modified gravity theories can 
be reproduced, at least for small deviations from GR, by 
a perturbative expansion around a background solution 
given by the standard TOV equations, which determine 
the structure of relativistic stellar models in GR [33, 34]. 

Instead of committing to one particular pet theory, 
in this paper we formulate a parametrized “post-TOV” 
framework for relativistic stars based on the well-known 
parametrized post-Newtonian (PPN) theory developed 
by Nordtvedt and Will [35, 36]; see e.g. [34, 37] for in¬ 
troductions to the formalism. The foundations of post- 
Newtonian (PN) theory for fluid configurations in GR 
were laid in classic work by Ghandrasekhar and collabo¬ 
rators [38, 39]. Various authors studied stellar structure 
using the PN approximation, both in GR [40-43] and in 
modified theories of gravity, such as scalar-tensor theory 
[44, 45]. To our knowledge, after some early work that 
will be discussed below [46-48] , the investigation of com¬ 
pact stars within the PPN approximation has remained 
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^2 = TTi± 7r2 ^ ± TTar^p ± 7r4 —, (2c) 
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Here r is the circumferential radius, m is the mass func¬ 
tion, p is the fluid pressure, p is the baryonic rest mass 
density, e is the total energy density, and H = {e — p)/p 
is the internal energy per unit baryonic mass. A “GR” 
subscript denotes the standard TOV equations in GR 
[cf. Eq. (7) below, where we appended a subscript “T” 
to the mass function for reasons that will become appar¬ 
ent later]; 6i,TTi (i = 1, ... ,4) and pi (i = 1,... ,5) are 
phenomenological post-TOV parameters. The GR limit 
of the formalism corresponds to setting all of these pa¬ 
rameters to zero, i.e. 5i,TTi,pi -T 0. 
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FIG. 2. 2PN-order post-TOV corrections on the mass-radius curves. We show the modification induced by different families of 
post-TOV terms on the general relativistic mass-radius curve, assuming the APR EOS. Left to right and top to bottom, the 
different panels show the effect of the pressure terms, proportional to (i = 1,... ,4), and of the mass terms, proportional to 
p,i (j = 1,... ,5). 


The dimensionless combinations Vi,A4i and 7^2, 2 

represent a parametrized departure from the GR stellar 
structure and are linear combinations of IPN- and 2PN- 
order terms, respectively. In particular, the coefficients 
Si attached to the IPN terms are simple algebraic combi¬ 
nations of the traditional PPN parameters: see Eqs. (35) 
and (36) below. As such, they are constrained to be very 
close to zero by existing Solar System and binary pulsar 
observations^: |(5i| ^ 1. This result translates to negli¬ 
gibly small IPN terms in Eq. (1): Vi,A4i 1. On the 


we 


other hand, tt^ and fit are presently unconstrained, and 
consequently 1 ^ 2 , M .2 should be viewed as describing the 
dominant (significant) departure from GR. 

Each of the two combinations 7^2 and Ad 2 involves no 
more than five dimensionless 2PN terms, but as we show 
in Section IIIB these terms are representative of five 
distinct “families” consisting of a large number of 2PN 
terms. Each family is defined by the property that all of 
its members lead to approximately self-similar changes in 


obtain the following upper limits: |<5i| < 6 X 10 ^,|52| $ 
7 X 10-3, l^al < 7 X 10-3, 1^41 < 10-®. 


^ Using the latest constraints on the PPN parameters [49] 
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the stellar mass-radius curves when included in 7 ^ 2 , 2 - 

In other words, as we verified by numerical calculations, 
we can account for several terms belonging to the same 
family by taking just one term from that family (either 
the dominant one or, when convenient, a much simpler 
subdominant one) and varying the corresponding post- 
TOV coefficient tt^ or fii. 

The qualitative effect of each of the 2PN-order post- 
TOV terms on the mass-radius relation is illustrated in 
Fig. 2. The values of the tt^ and /i^ coefficients in each 
panel of this figure were chosen for purely illustrative pur¬ 
poses, i.e., we chose these coefficients to be large enough 
that they can produce visible deviations on the scale 
of the plot. A first noteworthy feature is that pressure 
terms typically induce corrections that are about an or¬ 
der of magnitude smaller than mass terms^. This can be 
seen by the larger range of tt^’s needed to produce visible 
changes in the mass-radius curve (| 7 r 2 | < 4, IttsI < 100 
and | 7 r 4 | < 10 ) when compared to the corresponding cor¬ 
rections in the mass-function equation (|/i 2 | < 1 , l/ial < 1 
and 1/^41 < 1.5, respectively). Some terms (such as those 
proportional to 7 r 2 , TTa, 7 r 4 , ^3 and induce smooth 
rigid shifts of the mass-radius curve, similar to those 
that would be produced by a softening or stiffening of 
the nuclear EOS. Other terms (like those proportional to 
fii, p 2 and Pi) produce more peculiar features that are 
more or less localized in a finite range of central densi¬ 
ties. This is interesting, because (for example) it is plau¬ 
sible to conjecture that some combination of the pi and 
P 2 corrections may reproduce the qualitative features of 
a highly non-perturbative phenomenon such as sponta¬ 
neous scalarization, despite the intrinsically perturbative 
nature of our formalism. 

The punch line here is that each post-TOV correction 
is qualitatively different, so we can use the post-TOV for¬ 
malism as a toolbox to reproduce the mass-radius curves 
shown in Fig. 1 for various modified theories of grav¬ 
ity. More ambitiously, it would be interesting to address 
the inverse problem, i.e. to find out how the post-TOV 
parameters are related to the dominant corrections in¬ 
duced by each different theory. These issues are beyond 
the scope of this paper, but they are obviously crucial 
to relate our formalism to experiments, and we plan to 
address them in future work. 

The second main result of this paper has to do with 
the “completeness” of our post-TOV formalism, in the 
sense that the stellar structure Eqs. (1) - if we neglect 
the small terms - can be formally derived by a 

covariantly conserved perfect fluid stress energy tensor. 
That is: 

= 0, = (Ceff + p)u^u’' + pg^’', (3) 

where the effective, gravity-modified energy density is 

Ceff = e + pAf 2) (4) 


and the covariant derivative is compatible with the effec¬ 
tive post-TOV metric 

=diag[e''^’'\(l-2m(r)/r)“\r^,r^sin^6»], (5) 


with 


du 

dr 




( 6 ) 


Our phenomenological post-TOV formalism is expected 
to encompass a large number of alternative theories of 
gravity, but it is not completely general, and future ex¬ 
tensions may be possible or even desirable. As we stated 
earlier, theories which produce non-perturbative phase 
transitions in their stellar structure equations may not 
be accurately modeled. The formalism is also limited by 
the choice of acceptable 2PN terms out of all dimension- 
ally possible combinations, based on criteria that have 
bearing on the structure of the gravitational field equa¬ 
tions (see Section III B below). 


B. Plan of the paper 

The plan of the paper is as follows. In Section II we in¬ 
troduce the PPN formalism and review previous applica¬ 
tions to relativistic stars (in particular work by Wagoner 
and Malone [46] as well as Ciufolini and Ruffini [47]). In 
Section III we develop the post-TOV formalism to IPN 
order (where all parameters are already constrained to 
be very close to their GR values by Solar System and 
binary pulsar experiments), and then to 2PN order. We 
also show the equivalence between the 2PN post-TOV 
equations and GR with a gravity-modified EOS under 
a minimal set of reasonable assumptions. In Section IV 
we present some numerical results illustrating the rela¬ 
tive importance of the different post-TOV corrections. 
Some technical material is collected in three appendices. 
Appendix A gives details of the dimensional analysis ar¬ 
guments used to select the relevant set of 2PN post-TOV 
coefficients. In Appendix B we present a brief summary 
of the relativistic Lane-Emden theory, which plays an 
auxiliary role in the construction of our formalism. Fi¬ 
nally, Appendix C shows that certain integral potentials 
appearing at IPN order in the stellar structure equations 
(namely, the gravitational potential U, the internal en¬ 
ergy E and the gravitational potential energy O) can be 
approximated by linear combinations of non-integral po¬ 
tentials, so these integral potentials are “redundant” and 
can be discarded when building our post-TOV expansion. 


II. SETTING THE STAGE: STELLAR 
STRUCTURE WITHIN PPN THEORY 

A. The TOV equations 


^ A notable exception to this rule is the tti term, for reasons that 
will be explained in Section IV below. 


A convenient starting point for our analysis is the stan¬ 
dard general relativistic TOV equations, describing hy- 
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drostatic equilibrium in spherical symmetry [34]. These 
are given by the familiar formulas: 


dp\ _ (e + p) (tot + 47rr3p) 

r2 ' (1 - 2rm/ry^ ^ 


dr J, 


(7b) 


GR 


where p and e are the fluid’s pressure and energy den¬ 
sity, respectively, and tot is the mass function (the sub¬ 
script is used to distinguish this mass function from sim¬ 
ilar quantities appearing in PPN theory, see below). 

For later convenience we also write down the IPN- 
order expansion of these equations (for simplicity the 
subscript “GR” is omitted): 


J (l + n + ? + + 0 ( 2 PN), 

( 8 a) 

= 47rr^p(l -I- n). (8b) 

dr 

where we have introduced the baryonic rest-mass density 
p and the dimensionless internal energy per unit mass, 
n = (e — p)/p. It can be noticed that the mass function 
equation only contains IPN corrections to the Newtonian 
equations of hydrostatic equilibrium, while higher-order 
corrections appear in the pressure equation. 


B. The PPN stellar structure equations 


The PPN formalism [35, 36] was first employed for 
building static, spherically symmetric models of compact 
stars by Wagoner & Malone [46], and subsequently by 
Ciufolini & Ruffini [47]. This early work is briefly re¬ 
viewed here since it will provide the stepping stone to¬ 
wards formulating our post-TOV equations. 

A convenient starting point is the set of stellar struc¬ 
ture equations derived in [47] from the original Will- 
Nordtvedt PPN theory [35, 36]. These are [cf. Eqs. (14) 
of [47]]: 


dp 

dr 


em 

„2 


1 -b (5 -f 37 - 6^ -b C 2 ) — + - + Ca ? 

re m 


+ (7 + wi^ + i(n + 7-12^ + <.-2a)? 

TO 2 TO 


dm 

dr 


= 47rr^e, 


(9a) 

(9b) 


where we have adopted the standard notation for the 
nine PPN parameters, {/3, 7 , Ci, C 2 ) Ca, C 4 , ai, 02 , as}- In 
the GR limit /3 = 7 = 1 and (^i = ai = 0 {i = 1, ...,4) 
[49]. 

It should be pointed out that the basic parameters p, m 
(as well as the radial coordinate r) entering Eqs. (9) may 


not be the same as the corresponding ones in the TOV 
equations. This is a reflection of the “gauge” freedom in 
defining these parameters in a number of equivalent ways. 
Indeed, below we are going to exploit this freedom and 
obtain an “improved” set of PPN equations by a suitable 
redefinition of the mass function. On the other hand, fol¬ 
lowing [47] , we will stick to the same p and r throughout 
this analysis, implicitly assuming that they are the same 
variables as the ones in the TOV equations (7). 

The potentials and E appearing in Eq. (9a) obey 

— = —Airrpm, —— = 47 rr pH. (10) 
dr dr 

The more familiar Newtonian gravitational potential U, 
solution of = — 47 rp, is not featured in Eqs. (9) as a 
result of a change of radial coordinate and a redefinition 
of the mass function fh with respect to the original PPN 
theory parameters (see [47] for details). 

The stellar structure equations can be manipulated 
further by switching to a new mass function: 

ifi^ 

m(r) = fh + AE -b BQ + C -b DiAirr^p), (11) 

r 

where A, B, C, and D are free constants. As evident, fh 
and TO differ at IPN level. The constants A and B can 
be chosen so that the terms proportional to E and il in 
Eq. (9a) are eliminated. This is achieved for 

A = C3, 5= ^(11+7-12^ + C2-2C4). (12) 


The resulting “new” set of PPN stellar structure equa¬ 
tions is 


dp 

dr 


pm 

„2 


l + II-l- b(5-b 37 — 6(3 -b C 2 ~ C) — 

p r 


+ il + Ci- D) 47r 


r^p 


(13a) 


dm 

dr 


^ = 47rr^p 1-b (1-b Ca)!!-b 3£)-- 

[ p Air pr^ 

1 m 

-- (11 + 7 - 12/3 + C 2 - 2C4 -4C + 2 D) - 

2 r 


(13b) 


These expressions still contain the gauge freedom as¬ 
sociated with the definition of the mass function to in 
the form of the yet unspecified constants C and D. In 
particular, the Wagoner-Malone hydrostatic equilibrium 
equations [46] represent a special case of these expres¬ 
sions, and it is straightforward to see that they can be 
recovered for 


Z7 = 7 + C4, C= i(7 + 37-8/3 + C2). (14) 

Making this choice for the constants on the right-hand 
side of Eq. (11) leads to a new mass function, say to. 
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and to the following structure equations, which match 
Eqs. (6) and (7) of [46]: 


dp 

dr 

dm 

dr 


pm 

j,2 


^ ^ p m 

1 + n H-h a— 

p r 


(15a) 


p 


1 + (1 + Ca)!! + a -h 8(7 + C4)- 

r p 


b m? 

47r pr'^ ’ 


(15b) 


where a = (3 + 87 — 4/3 + C2)/2 and the constant b in the 
notation of [46] is our C, i.e. b = (7 + 87 — 8/3 + (,2) /2. 

A comparison between the two sets of PPN equations 
(9) and (15) discussed in this section reveals that the 
Wagoner-Malone equations are simpler, in the sense that 
they do not depend on the auxiliary potentials and E. 
This advantage, however, is partially offset by the more 
complicated expression for the mass function equation. If 
we compare the GR limit of the Wagoner-Malone equa¬ 
tions (15) against the IPN expansion of the TOV equa¬ 
tions, Eqs (8), we find that the two sets coincide provided 
we identify m = mx, i-C. 


m = mx 


O 

r 


+ Airr^p , 


(16) 


where the last equation follows by taking the GR limit of 
Eq. (11) in combination with Eqs. (12) and (14). Glearly, 
the fact that rh ^ tox in the GR limit is an unsatisfactory 
property of the Wagoner-Malone equations. 

It would be desirable to have a set of structure equa¬ 
tions that - unlike the set (9) - does not involve integral 
potentials, and such that - unlike the set (15) - the mass 
function is compatible with the GR limit. Fortunately, it 
is not too difficult to find a new set of PPN equations for 
which m = tox- In the following section we will propose 
an improved set of PPN stellar structure equations that 
satisfies these requirements. 


C. An improved set of PPN equations 


We can exploit the degree of freedom associated with 
the constants C, E in Eqs. (13) and produce a new set of 
PPN equations that exactly match the IPN TOV equa¬ 
tions in the GR limit with m = mx • It is easy to see that 
this can be achieved by making the trivial choice 


Note that the constants A and B are still given by 
Eqs. (12). The resulting PPN equations are 


dp 

dr 


dm 

dr 


pm 

r.2 


T) Tfl 

1 -I- n -I- - -I- (5 -I- 87 — 6/3 -I- C 2 ) — 

r 


-b (7 -f C4)47r 


P 

Q 

r p 
m 


(18a) 


47rr^p[l -b (1 -b Ca)!! 

1 777 

— - (11 -b 7 — 12/3 -b C2 — 2^4) — 

2 r 


(18b) 


As advertised, in the GR limit these equations reduce to 
Eqs. (8) with m = rn-x- The same equations will be used 
in Section III below in the construction of the desired 
post-TOV equations. 


D. The physical interpretation of the mass function 

Within the framework of PPN theory, inertial mass 
and active/passive gravitational mass are, in general, dis¬ 
tinct notions. In the context of compact stars, expres¬ 
sions for all three kinds of mass are given in [47]: 

= m{R) + (^^ + ^7 - 10/3 + ^C2) m), (19) 

Ma = Mill + (^4/3 — 7 — 3 — -tta — — 2^2^ VL{R) 

+ Q:iE(R) - Qaa - 3C4 + Ci) ^, (20) 

Mp = Min -b (^4/3 - 7 - 3 - ai -b ^02 - ^Ci - ^(2^ 

X fl(i?), (21) 

where R is the stellar radius associated with the mass 
function fh(r) - i.e. with the set of equations (9) - and 

nR 

P = 4tt drr^p. (22) 

Jo 

is the volume-integrated pressure. 

In GR the three masses are of course identical, Min = 
Ma = Mp. As argued in [47], any theory conserving the 
four-momentum of an isolated system should incorporate 
the equality of the two gravitational masses, i.e. Ma = 
Mp. If adopted, this equality leads to following three 
algebraic relations for the PPN parameters: 

Ca = 0, (23) 

Cl — 3C4 + 2*^3 = 0, (24) 

Cl + 3ai — 2 q!2 ~ 5C2 — 2®3 = 0. (25) 

We can subsequently write for the common gravitational 
mass: 


C = D = 0. 


(17) 


Mg = Ma = Mp = m{R) + ra(i?), 


(26) 
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with 


A. Post-TOV equations: IPN order 


i^=^(ll + 7-12/3-a3 + C2-^Ci)- (27) 

For our new PPN equations with C = D = 0 the mass 
equality Ma = Mp implies 

m{r) = fh{r) + ^ (11 + 7 - 12/3 + C2 - 2(4) ^2(r). (28) 

Then with the help of Eq. (24) it is easy to see that 

Mg = m{R) + ^(.R) = m{R). (29) 

If R is the stellar radius associated with our PPN equa¬ 
tions (18), the difference SR = R — R is a IPN-order 
quantity. We can then approximately write 

m{R)^m{R)-^iR)SR. (30) 

dr 

However, Eq. (18b) implies that dm/dr{R) = 0 if p{R) = 
0 at the stellar surface. This is indeed the case for a 
realistic EOS. Therefore, we have shown that at IPN 
precision the mass of the system is given by 

Mg = m(i?). (31) 


The recipe for formulating leading-order post-TOV 
equations is rather simple: from a suitable set of PPN 
hydrostatic equilibrium equations we isolate the purely 
non-GR pieces. These IPN terms are subsequently added 
“by hand” to the full general relativistic TOV equations, 
hence producing a set of parametrized post-Einsteinian 
equations. It should be pointed out that this proce¬ 
dure can only be applied at the level of IPN correc¬ 
tions. Higher-order corrections should by sought by other 
means, such as dimensional analysis (see Section HIB). 

In principle, either set of equations, (9) [47] or (15) 
[46], could have been used. However, our improved PPN 
equations (18) seem to be best suited for this task. 

Considering Eqs. (18), we first isolate the terms that 
represent a genuine deviation from GR. These are the 
terms in the second line in the following equations: 


dp 


pm 


p 2m 


■3 

r p 


l + +-+47r — 


pm 

^2 


p r m 

Tfi T D 

(3 -I- 37 - 6/3 + C2)-1- (7 - 1 + C4)47r- 

r m 


dm , 2 

—— = 47rr^p (1 -I- H) 
dr 


+ Anr^p 


1 


m 


C3n--(ll + 7-12/3 + C2-2C4) 

2 r 


(32a) 


.(32b) 


This elegant result is one more attractive property of the 
new PPN equations. 


III. THE POST-TOV FORMALISM 


The logic underpinning the formalism we are seeking 
is that of parametrizing the deviation of the stellar struc¬ 
ture equations from their GR counterparts, thus produc¬ 
ing a set of post-TOV equations. As already pointed out 
in the introduction, the post-TOV formalism is merely a 
useful parametrized framework rather than the product 
of a specific, self-consistent modified gravity theory (in 
the spirit of PPN theory). In this sense our formalism is 
akin to the existing “quasi-Kerr” or “bumpy” Kerr met¬ 
rics, designed to study deviations from the Kerr space- 
time in GR (see e.g. [51-53]). 

By design the post-TOV formalism should be a more 
powerful tool for building relativistic stars than the PPN 
framework; after all, the latter is based on a PN approx¬ 
imation of strong gravity. However (as will become clear 
from the analysis of this section), our formalism has its 
own limitations, the most important one being the fact 
that the deviations from GR are introduced in the form 
of PN corrections. This could mean that the structure of 
compact stars with a high degree of departure from GR 
may not be accurately captured by the formalism. 


The second step consists of adding the non-GR terms to 
the TOV equations (7). We obtain (recall that m = mx) 

m + Anr^p \ pm 
1 — 2m/r J 


dp ^ je + p) 
dr r^ 


Si- 




Q 

r p 
m 


dm , 9 

—j— = Airr 
dr 



+ <34n^ , 


(33) 

(34) 


where we have introduced the constant post-TOV pa¬ 
rameters: 


= 3(1-b7) - 6/3-f C2, (52 = 7-1 + C4, (35) 

^3 = -^ (11 + 7 - 12/3 + C2 - 2C4), <34 = C3- (36) 

As expected, (5^ = 0 in the limit of GR. 

The above equations can be written in a more compact 
form: 


dp 

fdp\ 

dr 

\drj 

dm 

/ dm^ 

dr 

V dr , 


GR 

^GR 


pm 

„2 



-I- i52 47r 



-I- 47rr^p '^411) • 


(37a) 

(37b) 


These expressions represent our main result for 
the leading-order post-TOV stellar structure equations. 
They describe the IPN-level corrections produced by an 
arbitrary deviation from GR that is compatible with 










PPN theory. In other words, Eqs. (37) encapsulate the 
stellar structure physics (at this order) for any member 
of the PPN family of gravity theories. 

We could in principle introduce other IPN order terms, 
in the spirit of the general parametrized framework of 
deviating from GR that we have described in the begin¬ 
ning of this section. But the introduction of such terms 
would correspond to either redefinitions of coordinates 
and/or the mass function at the IPN level, as we have 
already seen, or deviations from special relativity, which 
we would prefer not to include. 

Unfortunately, it turns out that Eqs. (37) are of lim¬ 
ited practical value. As discussed in the executive sum¬ 
mary, the modern limits on the PPN parameters suggest 
that these corrections are very close to their GR values, 
because /3 ,7 « 1 and OiXi 1 ; making all the 6i param¬ 
eters very small. We should not therefore expect any no¬ 
table deviation from GR at the level of the leading-order 
post-TOV equations. We verified this claim by explicit 
calculations of neutron star stellar models with different 
EOSs. 

Any significant deviations from compact stars in GR 
have to be sought at 2PN order and beyond, where the 
existing observational limits leave much room for the 
practitioner of alternative theories of gravity. This calls 
for the formulation of a higher-order set of post-TOV 
equations, a task to which we now turn. 


B. Post-TOV equations: 2PN order 


In this section we shall formulate post-TOV equations 
with 2PN-accurate correction terms. Unlike the calcu¬ 
lation of the preceding section, we now have to build 
these equations “from scratch”, given that the general 
PPN theory has not yet been extended to 2PN order. 
Inevitably, the procedure for building the various 2PN 
terms will turn out to be somewhat more complicated 
than that of the preceding section, heavily relying on di¬ 
mensional analysis for constructing these terms out of 
the available fluid parameters. Moreover, at 2PN order 
we also need to consider terms that involve the integral 
potentials U, E and U (recall that these were eliminated 
at IPN order by a suitable redefinition of the mass func¬ 
tion). However, as shown numerically and via analytical 
arguments in Appendix G, the integral potentials can be 
approximated to a high precision, and for a variety of 
EOSs, by simple linear combinations of the non-integral 
PN terms. As a result, they do not have to be considered 
separately in the post-TOV expansion. 

To begin with, we can get an idea of the form of some 
of the 2PN terms we are looking for by expanding the 
TOV equations (7) to that order. Let us first consider 


the pressure equation (7a): 


dp 

dr 


pm 

J.2 





2m , r^p 

1 -\ -h 47r- 

r m 


Am? „ 9 

-I- Y~ + 87rr^p 


-ke)(3PN). 


( 38 ) 


As anticipated, all IPN terms appearing here are also 
present in our PPN equation (I 8 a). The produced 2PN 
corrections are proportional to the following combina¬ 
tions: 


-, H —, r^p, 


n 


“^9 

r p r p 


rp 


pm 


(39) 


Additional 2PN terms that do not appear in the TOV 
equations can be constructed by forming products of the 
available IPN terms. The largest set of IPN terms can 
be found in the general PPN equations (13): 


IPN : 


n, 


p m r p 


m 



(40) 


We can observe that all terms, except the last one, also 
appear in our final PPN equations (18). From these we 
can reproduce the set (39) as well as the additional 2PN 
terms: 


m^ ’ pr'^ ’ 


m 


pr° 






p 


m 


m?p 


p2 ’ p2^S. ’ p2^i ■ 


(41) 


We have set apart the last three (underbraced) terms of 
this set because, as a result of their ^ 1 /p^ scaling, these 
terms will be discarded. In fact, the same fate will be 
shared by any term ~ p^ with /3 < —2. 

There are various reasons why we believe that this se¬ 
lection rule should be imposed. In our opinion these rea¬ 
sons are quite convincing, but they fall short of consti¬ 
tuting a watertight argument: in all fairness, if we had a 
single truly compelling reason, we would not need more 
than one. 

The first line of reasoning to exclude the presence of 
negative powers of p (and of the other fluid parameters) 
in the PN terms is based on the regularity of these terms 
at the stellar surface, where p, p, H —)• 0 for any realistic 
EOS. A PN term like the second one in the underbraced 
group of the set (41) will lead to a term diverging as 
~ 1 /p at the stellar surface in the stellar structure equa¬ 
tions, and therefore it is not an acceptable PN correction. 
Although this surface regularity argument is powerful, it 
obviously works only for terms that do not scale with 
positive powers of p or H. 

The second - heuristic - argument applies to gravity 
theories with the following (symbolic) structure: 


{geometry} = 8nT^^, 


(42) 


= 0 


dp 

dr 


(e+p){ geometry}. 


(43) 
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where “geometry” stands for combinations of the met¬ 
ric and its derivatives, and the last equation assumes a 
perfect fluid stress-energy tensor. The stress-energy ten¬ 
sor and the right-hand side of Eq. (43) feature e + p = 
p{\ -\-IV + p/p) and p linearly. It can then be argued that 
the solution of the field equations for the metric and its 
derivatives will display a 

{geometry} - (e -h rp)" ~ p” (^1 -b H -H (44) 

dependence with respect to the fluid variables (where t 
and n are 0{1) numbers). Such a solution should lead to 
pressure-dependent PN terms of the form: 

PN term ~ (r^p)"“^ , k = n,n—l,... (45) 

where one p factor has been removed and absorbed in the 
Newtonian prefactor of the structure equations, while at 
the same time the factor has been added in order to 
produce a dimensionless quantity. A key observation is 
that the form (45) assumes a theory that does not de¬ 
pend on dimensional coupling constants. Now, accord¬ 
ing to (45) the highest negative power of p corresponds 
to k = n, which means that the scaling with respect to 
the density should be: 

PN term ^ p^ ^ /3 > — 1 (46) 

Based on these arguments, we deem acceptable those PN 
terms which scale with p as in (46). This choice is also 
consistent with the previous PPN formulas, see Eqs. (13). 
A similar argument can be used to exclude terms with 
negative powers of p and 11^. 

Equation (39) and the top row of Eq. (41) represent 
a large set of 2PN terms emerging from the expansion 
of the TOV equation and from products of the various 
known IPN terms. This set is large but not necessar¬ 
ily complete. Inevitably, a systematic approach to the 
problem of “guessing” 2PN terms should involve dimen¬ 
sional analysis. To improve readability we relegate our 
dimensional analysis considerations to Appendix A, and 
here we only quote the main result. The most general 
form for 2PN order terms is given by the dimensionless 
combination: 

A2 ~ n^r^prir^pf (™) , (47) 


^ A related argument for excluding high powers of 1/p is the fol¬ 
lowing. By virtue of the field equations, the Ricci scalar is usu¬ 
ally proportional to the energy density of matter (at least in 
the Newtonian limit, if the modified theory reproduces GR in 
the weak field regime): R ~ p. If inverse powers of p are pro¬ 
duced by gravity modifications, they should therefore originate 
from terms ~ in the action of the theory. These terms are 

usually associated with ghosts or instabilities [54], and therefore 
their presence is problematic. 


where a, /3, ^ are integers with 


/?>-!, 


(48) 


while different bounds on 9 and a apply to the two hy¬ 
drostatic equilibrium equations: 


0<6»<2, 0<a<2-6», (49) 

dr 

UjTTL 

— : 0<6'<3, 0<q;<3-6». (50) 

dr 

The lower bounds on the three parameters a, /S, 9 
are dictated by the same considerations discussed below 
Eq. (41), namely, regularity at the surface and consis¬ 
tency with the fact that gravitational field equations of 
the general form (43) are unlikely to generate negative 
powers higher than 1/p. The upper bounds on a and 
9 are imposed by the regularity at r = 0 of the stellar 
structure terms arising from A2 (see Appendix A). 

From the general expression (47) we can reproduce all 
previous IPN and 2PN terms and generate an infinite 
number of new ones. This possibility could have been a 
fatal blow to our post-TOV program. Fortunately, the 
day is saved by the fact that the magnitude of A2 decays 
rapidly throughout the star as /3 increases. This trend 
is clearly visible in the numerical results shown in Fig. 3 
(see discussion below). 

For all practical purposes these results imply that the 
first few members of the /3 = —1,0,1,... sequence are 
sufficient to construct accurate post-TOV expansions. A 
sample set of such dominant 2PN terms is: 


r5p 


mp 2 r^p^ 
y, rpm, r p, - 


2PN : rpm, — 


rp 

,2 


pm 


rSp2 


7 10 “^ 

r P m‘ m 2 rrP 

—T —n^, n—, Hr p, n- 

r^p r 


pm'^ 


P 

^P ^2 "1 ^2 n2 

11 „ , 11 „ , 11 , 11 , 

pr'^ mp 


n—,n 

m pm 
^2 r^P 




9 1 

r^p 


m 


(51) 


This set is markedly larger than the previous sets (39) 
and (41) (whose acceptable terms form a subset of the 
new set), but a complete post-TOV formalism would have 
to include all (or almost all) of these terms, with twice 
the number of free coefficient in the dp/dr and dm/dr 
equations. Fortunately, as it turns out, the same job 
can be done with a much smaller subset of 2PN terms. 
This is possible because the various 2PN terms can be 
divided into five “families ”, each family comprising terms 
with similar profiles. When incorporated in the post- 
TOV equations, terms belonging to a given family lead 
to self-similar modifications in the mass-radius curves for 
a given EOS. 

Insight into the behavior of the h. 2 {a, 13,9) terms can 
be gained by direct numerical calculations of their radial 
profiles in relativistic stars. We carried out such calcula¬ 
tions for a variety of realistic EOSs as well as relativistic 
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r/B. 






FIG. 3. The radial profile of the A. 2 -function. We exhibit the behavior of pA. 2 , where p = p/ pc and A 2 is given in Eq. (47), 
for a stellar model using the APR EOS, with ede? = 0.86 x 10^® g/cm®, M = I.SIMq and R = 12.3 km. The curves are 
labelled according to the respective values of (a,/3,6). From the top row to the bottom row the index 9 takes on the values 
(0,1, 2, 3), respectively. Despite the multitude of possible dimensionally correct 2PN terms, their self-similarity - which is clear 
when we compare terms along the bottom-left to top-right diagonals in this “grid” of plots - allows us to group them into a 
relatively small number of families (see text for details). The contributions plotted in three panels at the bottom right of the 
grid (marked as “Excluded”) would lead to divergences in the hydrostatic equilibrium equations, and therefore they can be 
discarded as unphysical. 


polytropes, and for different choices of central density, 
verifying that all cases lead to very similar results, as 
discussed below. More specifically, we considered EOS 
A [55], EPS [56], SLy4 [57, 58] and N [59] in increasing 
order of stiffness, as well as relativistic polytropes with 
indices n = 0.4, 0.6, and 1.0: see Appendix B, and in 
particular Eq. (B7). Note that the polytropic models are 
parametrized by A = Pc/cc instead of Cc alone (the sub¬ 
script “c” indicates a quantity evaluated at the center), 
but this is equivalent to the central density parametriza- 


tion. The polytropic models are also invariant with re¬ 
spect to the scale factor this can be adjusted to 

generate polytropic models of (say) the same mass (for a 
given A) as that of a specific tabulated-EOS model. 

Rather than computing A2 itself, from a phenomeno¬ 
logical point of view it makes more sense to consider 
the combination pK 2 - The reason is that this combi¬ 
nation appears in both the pressure and mass equa¬ 
tions, and furthermore it has the desirable feature of 
being regular at the surface for (3 = —1. More specif- 
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FIG. 4. The family-representative 2PN terms. Here we show the selected representative terms from each of the families depicted 
in Fig. 3, as listed in Eq. (54), for three different EOSs: FPS (left panel), APR (middle panel) and an n = 0.6 polytrope (right 
panel). Each term illustrates the qualitative behavior of each family of possible 2PN contributions to the structure equations. 
The high degree of invariance of the A 2 -profiles with respect to the EOS is evident in this figure. The GR background stellar 
models utilized in the figure have the following bulk properties: tc = 0.861 x 10^® g/cm^ (A = Pc/cc = 0.165), M = 1.51 M© 
and R — 12.3 km (left panel); Cc = 1.450 x 10^® g/cm® (A = 0.198), M = 1.50 M© and R = 10.7 km (center panel); A = 0.165, 
M = 1.50Mq and R — 11.75 km (right panel). 


ically, in Fig. 3 we plot the dimensionless combination 
pA 2 (a, P,0), where p = p/pc- Our sample neutron star 
model was built using the APR EOS with central energy 
density tc/e? = 0.86 x 10^® g/cm^, corresponding to mass 
M = I.5IM0 and radius R — 12.3 km, but we have veri¬ 
fied that our qualitative conclusions remain the same for 
different models and different EOSs. 

Eigure 3 reveals two key trends: (i) the clear /3-ordering 
of the pA 2 profiles, with /3 = — 1 always associated with 
the dominant term for fixed a and 9, and (ii) the re¬ 
markable similarity in the shape of the profiles of terms 
with dissimilar (a, /3, 9) triads along the bottom-left to 
top-right diagonals in the “grid” of Fig. 3. This prop¬ 
erty defines distinct families of 2PN terms and implies 
that the terms of each family cause self-similar changes 
in the mass-radius curves of the various post-TOV stellar 
models. 

We have identified five 2PN families (labeled 
“F1”,...,“F5” in the various panels of Fig. 3, and de¬ 
scribed in more detail in Table I): 

(i) FI: This is a single-member family comprising only 
the pA 2 { 0 , —1, 0) term in the top-left panel, which is zero 
at r = 0 but finite a.t r = R. 

(ii) F2: The members of this family vanish at r = 0 
and r = R, and have a peak near the surface. These are 
the pA2(0,/3,0) terms with /? > 0 in the top-left panel. 

(hi) F3: These terms also vanish at both r = 0 and 
r = R, but display an approximately flat profile inside 
the star. They correspond to pA2(l,/3,0) (top-middle 
panel) and pA2(0,/3,1) (top-right panel) for /3 > —1. 

(iv) F4: This family comprises terms that are fi¬ 
nite at r = 0 but zero at r = i?. These are the 
pA2(2, /3, 0) (bottom-left panel) and pA2(l, /3,1) (bottom- 
middle panel) terms with /3 > — 1. 

(v) F5: These terms by themselves diverge at r = 0 
and vanish at r = ii, but they become well-behaved when 
inserted into the stellar mass-function equation, where 


they are multiplied by the factor r^: cf. Eq. (A24). 
These terms correspond to pA2(2, /3,1), and from the con¬ 
straints (49) and (50) we conclude that the members of 
this family can only appear in the mass equation. 

There is an intuitive way to explain the existence of 
the above families. As an example we consider E3, where 
the seemingly unrelated terms A2(l,/3,0) and A2(0,/3,1) 
yield similar profiles. Consider 

A2(0,/3,l)~r-i+3/5^. (52) 

By means of the approximations m ~ pr^, If ~ p/p (the 
latter approximation is motivated by the exact thermo¬ 
dynamical relation fl = np/p for relativistic polytropes 
with index n, see Appendix B) we find 

A2(0, /3,1) ~ ^ A2(1, /3,0). (53) 

\m/ 

Similarly we can show that A2(2,/3, 0) ~ A2(l,,0,1) for 
the F4 family. The argument can be generalized to show 
that terms along the diagonals of Fig. 3 are equivalent. 

Table I summarizes the taxonomy of the most impor¬ 
tant terms of each family according to the above criteria. 
The impact of each of these terms as a post-TOV correc¬ 
tion has been tested for a variety of EOSs. The results 
reveal that the members of a given family lead to self¬ 
similar modifications to the stellar mass-radius curves. 
A sample of these numerical results is shown in Fig. 2, 
which is further discussed in Section IV below. 

This remarkable self-similarity property means that we 
can simply select one term from each family and emulate 
the effect of all significant 2PN terms of the same family 
by simply varying the post-TOV coefficient associated 
with the selected term. 

In doing so, it is reasonable to choose the simplest 
terms as family representatives. For families F2 and F4 
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EOS APR 



Relativistic Polytrope (n =0.6) 




r/R 


FIG. 5. The family-representative terms in the structure equations. This figure illustrates the behavior of each of the family- 
representative 2PN terms [see Eq. (54)] multiplied by the Newtonian prefactors in the post-TOV equations. The stellar 
parameters are identical to the ones used in Fig. 4. Left panel: the combination {pm/r'^) A 2 {a, ft, 9) appearing in the pressure 
equation. Right panel: the combination pr^ A 2 {a, P, 6) appearing in the mass equation. The top panels correspond to EOS 
APR; the bottom panels correspond to a relativistic polytrope with pol 3 rtropic index n = 0.6. The divergence at the origin of 
the F5 term justifies its exclusion from the pressure equation. 


the simplest terms also happen to be the dominant ones 
(i.e., the ones with the largest ph- 2 ), while for F3 and 
F5 they are the first subdominant ones. The case of the 
single-member family FI is trivial. The five terms we 
select based on this reasoning are: 

TV? 2 P s ^ 

F-representatives : r p, 11-, 11 —. (54) 

r°p pm 

The phenomenologically relevant radial profiles of 
pA 2 {a, /3, 0) produced by these terms are shown in Fig. 4 
for three choices of EOS: FPS, APR and an n = 0.6 poly¬ 
trope. The most striking feature of this figure is the close 
resemblance of the A2 profiles of identical (a, /?, 6) tri¬ 
ads for different EOSs, which lends support to the EOS- 
independence of our selection of post-TOV terms. 

The family-representative terms (54) are again shown 
in Fig. 5, where we plot the combinations that appear 
in the dp/dr and dm/dr equations, i.e. mpA 2 /r^ and 
r^pA2, respectively (in the latter term we have omitted a 


trivial prefactor of 47r). We consider two different EOSs: 
APR and an n = 0.6 polytrope. All terms displayed are 
regular at both r = 0 and r = R with the exception of 
the F5 term in the dp/dr equation, which is divergent at 
r = 0 and must be excluded. Once again, the variations 
in the radial profiles due to considering different EOSs 
are extremely mild. 


We have thus obtained a minimum set of representa¬ 
tive 2PN terms, listed in Eq. (54), which in reality en¬ 
compasses a much larger set, like the one obtained from 
the combination of Eqs. (41) and (51), as well as terms 
that involve the integral potentials. 


After this admittedly tedious procedure we can fi¬ 
nally assemble our 2PN-order post-TOV equations for 
the pressure and the mass. These are (omitting the neg- 
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Family 

2PN term 

{a,p,e) 

Dominant / 
Chosen? 

El 

rA/{Ap) 

(o,-i,o) 

D/C 

F2 

{m/rf 

(0,0,0) 

D/C 

F2 

rmp 

(0,1,0) 

- 

F3 

mpl{rp) 

(l.-l,0) 

- 

F3 

r^p 

(1,0,0) 

- 

F3 

nm^/(r^p) 

(0,-l,l) 

D/C 

F3 

Ilm/r 

(0,0,1) 

- 

F3 

r^IIp 

(0,1,1) 

- 

F4 

r^p^/{pm) 

(2,-1,0) 

- 

F4 

r®p^/(m^) 

(2,0,0) 

- 

F4 

Up/p 

(l.-l.l) 

C 

F4 

Ilr^p/m 

(1,0,1) 

- 

F4 

H^m/(r^ p) 

(0,-1, 2) 

D 

F4 


(0,0,2) 

- 

F5 

nr'^p^/(pm^) 

(2,-1,1) 

- 

F5 


(2,0,1) 

- 

F5 

Yl^rp/mp 

(1,-1,2) 

- 

F5 

Tl^r^p/m^ 

(1,0,2) 

- 

F5 

nV(rV) 

(0,-1, 3) 

D 

F5 

n^r/m 

(0,0,3) 

C 


TABLE I. Taxonomy of the dominant 2PN terms. The self¬ 
similarity between the radial profiles of the various 2PN terms 
listed in Eq. (51) (and illustrated in Fig. 3) allows us to group 
them into five distinct families. This table spells out the ex¬ 
plicit form of the various terms, and indicates which term in 
each family is dominant (D) according to our numerical cal¬ 
culations, and which one was chosen (C) as representative of 
each family. 


C. Completing the formalism: the post-TOV 
metric and stress-energy tensor 

So far, our post-TOV formalism comprises no more 
than a pair of stellar structure equations, Eqs. (55a) and 
(55b), which can be used for the description of static and 
spherically symmetric compact stars. In this section we 
show that there is more to the formalism than meets the 
eye: to a high precision it is a “complete” toolkit, in the 
sense that (i) it can be reformulated in terms of a spher¬ 
ically symmetric metric and a perfect fluid stress- 
energy tensor T^’’', and (ii) these two structures are re¬ 
lated through the covariant conservation law = 0 

(where is the metric-compatible covariant derivative), 
hence respecting the equivalence principle. Remarkably, 
it also turns out that the metric and matter degrees of 
freedom can be related as in GR, which implies that the 
post-TOV formalism is equivalent to stellar structure in 
GR with a gravity-modified EOS for matter and an effec¬ 
tive spacetime geometry. 

In order to establish the above statements we begin 
with the following general result. Assume the static 
spherically symmetric metric 

ds^ = g^vdx^dx'" 

= -e<^)de +(^1- ' dr^ + r^dn^ (56) 

and a perfect-fluid stress-energy tensor (with energy den¬ 
sity £ and pressure V) 

= {£+ V)u>^u^+Vg^‘'. (57) 

Eor a static spherical fluid ball, the energy-momentum 
conservation equation 


ligibly small IPN corrections): 


= 0 (58) 

leads to 


dp ^ ( diA 
dr U^/gr 

+ TTar^p -b 7r4n- ) , (55a) 

P J 


pm 

-.2 


TTl 


j,5p 


■ 7^2- 


dm 

dr 


/ dm\ 
dr ) 


GR 


m 


m 


-bdTrr p[ + P2^ 


j.5p 


+ Tsr p -\- pilV —h 

p m 


(55b) 


where, as anticipated in the executive summary, tt^ (i = 
1,..., 4) and pi {i = 1,..., 5) are free parameters con¬ 
trolling the size of the corresponding departure from GR. 


^ = -(^£ + V)TU = -\{£ + V)^^^. (59) 

As long as we consider theories respecting (58) with a 
metric-compatible covariant derivative, this result is in¬ 
dependent of the gravitational field equations. 

Eor the mass function A4{r) we can always write a 
relation of the form 

= Airr'^S [1 -b Z{r)] , (60) 

dr 

where Z(r) is a theory-dependent function. Einstein’s 
theory is recovered by setting Z = 0, as required by the 
field equations of GR. 

To establish the properties described at the beginning 
of this section we will show that we can successfully map 
our post-TOV equations onto Eqs. (59) and (60) (with 
Z = 0). 
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The full post-TOV equations, Eqs. (55), can be written 
in the form: 


dp 

dr 


i^+P)T.f N P'm- 


1 + n + f^ 

p 


dm 

dr 


47rr^e + 47rr^p [ Af i + Af2 ] > 


Vi +V 2 , 
(61) 


(62) 


where Vi , V 2 have been dehned in Eqs. (2) , 


■p2 = 7^2 ~ 



(63) 


with 


dv 

dr 


^ {1 - M 2 )^ + rn(Vi+V 2 

^[{i-M2)r + mr2]. 


MiVi) 


(72) 


(73) 


From this result we can see that r represents the circum¬ 
ferential radius of the r = constant spheres and therefore 
the post-TOV radius R (where p{R) — 0) coincides with 
the circumferential radius of the star. 

Finally, the effective post-TOV stress-energy tensor is 


r'^" = (eeff+p)u'^u"+W^", (74) 


is a 2PN-order term, andr(r) = (m-|-47rr^p)/(l —2TO/r). 

Based on these expressions, we can define the effective 
energy density 

Eeff = e + p{Mi+ M 2 ), (64) 


and it is covariantly conserved with respect to the metric 
(71). 

These expressions clearly demonstrate that our post- 
TOV formalism is completely equivalent to GR with an 
effective EOS: 


which implies 


dm 

dr 


dTrr^Eeff- 


Using (64) in the pressure equation we have 


^ - [ieS + P - piMi +M 2 )] 4 

dr 

TfL r ~ 

- 2 (e + p)^^! + (Ceff + p)7^2 , 


(65) 


( 66 ) 


where we have used the fact that in any 2PN term we can 
replace the factor p by CeS -f p, at the cost of introducing 
3PN terms. Using (64) once more in the last term, and 
after some rearrangement, we obtain, 


^ _ (eeff -I- p) 
dr r^ 


+ 



(1 — Af 2)r -I- m{Vi + V 2 


MiVi) 

(67) 


Given that Af 1 <C 1, the last term can be safely omitted 
and we are left with 


dp 

dr 


(Ceff + P) 


(1 - Af 2 )r -b miVi +V 2 - MiVi) 

( 68 ) 


(Ceff + P) 


[(l-Af2)r-bmiP2] 


(69) 


which is of the form (59) . Note that in this and the follow¬ 
ing expressions the small Af 1 ,terms can be omitted. 
The resulting mapping is: 


V = p, M = m, S = Ces. (70) 

It follows that the effective post-TOV metric is 
Pf,,, = diag[(1 - 2m(r)/r)“\sin^ 6>], (71) 


p(e) ^ p(eeff), (75) 

eeff~e + pAf2- (76) 

As is evident from this last expression, Ceff represents a 
gravity-shifted parameter with respect to the physical en¬ 
ergy density e. This result highlights a key characteristic 
of compact relativistic stars when studied in the context 
of alternative theories of gravity, namely, the intrinsic de¬ 
generacy between the physics of the matter and gravity 
sectors. 

Whether the above effective description (and in par¬ 
ticular its effective geometry part) can give observables 
that have a correspondence to observables of an underly¬ 
ing theory or not depends on the nature of that theory. 
As long as the underlying theory admits a PN expansion, 
the physical description that arises from the effective for¬ 
malism should match that of the physical theory. This 
non-trivial issue will be further discussed elsewhere [60]. 


IV. NUMERICAL RESULTS 

In this section we provide a more detailed discussion 
of our numerical techniques and results, focusing on the 
mass-radius curves produced by the integration of the 
post-TOV equations (55) [or equivalently Eqs. (1)]. 

First, let us briefly summarize the integration proce¬ 
dure we have followed in this paper. We have carried 
out two kinds of computations: (i) “background” models 
- these involve the integration of the general relativistic 
TOV equations - with the purpose of studying the ra¬ 
dial prohles of the post-TOV correction terms, (ii) the 
integration of the full post-TOV equations, typically in¬ 
cluding the representative term of a single 2PN family. 

The post-TOV structure equations (la) and (lb) are 
integrated simultaneously starting at the origin r = 0, for 
fixed values of the coefficients tTj, pi, and for a range of 
central energy density values. The chosen central energy 
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FIG. 6. Self-similarity in mass-radius curves - I. Numerical integrations show that 2PN terms belonging to the same family 
result in self-similar deviations from GR in the mass-radius relation. This figure illustrates this remarkable property for pressure 
terms (top row) and mass terms (bottom row) belonging to families F2, F3 and F4 (from left to right). In each panel, the 
solid line corresponds to GR; the long-dashed (online: blue) line to a positive-sign correction due to the chosen term in each 
family; the short-dashed (online: red) line to a negative-sign correction due to the chosen term in each family. The various 
symbols show that nearly identical corrections can be produced using different terms belonging to the same family, as long as 
we appropriately rescale their post-TOV coefficients. 



R [km] 

FIG. 7. Self-similarity in mass-radius curves - 11. Same as in 
Fig. 6, but for the F5 family, which only admits post-TOV 
corrections with ^5 < 0 (see text). 


density Cc fixes the central pressure Pc = p(ec), the cen¬ 
tral mass density pc = TObnb(ec) and the central internal 
energy lie = (cc - Pc)IPc, where mb = 1.66 x g is 

the baryonic mass and rib is the baryon number density. 
In general, p{e) and nb(e) are computed using tabulated 


EOS data. Once the initial conditions have been speci¬ 
fied, Eqs. (la) and (lb) are integrated outward up to the 
stellar radius R, where p{R) — 0. The gravitational mass 
is obtained as M = m{R). 

The integration procedure for realistic EOS back¬ 
ground models is virtually the same as the one just de¬ 
scribed. We have also employed a number of polytropic 
background models; for these the integration procedure 
is slightly different (see Appendix B for details), and it 
is based on the simpler Lane-Emden formulation, where 
the pressure is replaced by the density p in the structure 
equations and the stellar model is parameterized by the 
ratio A = Pejsc rather than Cc alone (this formulation is 
of course equivalent to the one using tabulated EOSs). 
The added advantage of this approach is its scale invari¬ 
ance with respect to the polytropic constant K. This 
means that K can be freely adjusted to generate a model 
with (say) a specific mass M. This scaling procedure also 
fixes the radius R. 

The main installment of our mass-radius results has al¬ 
ready been presented in Fig. 2 of the executive summary 
(Section I A). As discussed there, the various post-TOV 
correction terms, representing the five 2PN families of 
Section IIIB, cause qualitatively different modifications 
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FIG. 8. Fractional deviations induced by the post-TOV parameters on the stellar mass and radius. Here we illustrate the 
fractional changes cansed by the post-TOV parameters in nentron star masses and radii. For a fixed central energy density 
and EOS APR, we plot the relative deviations from GR in mass and radius that result from varying the post-TOV parameters 
within the range indicated in the legends. Top row: e^jc? = 8.61 x 10^^ g/cm®, Mgr = 1.51 Mq and Rgr = 12.3 km. Bottom 
row: tc/c? = 1.20 x lO'^® g/cm®, Mgr = 2.04 Mq and 7?gr = 11-9 km. Left panels: Effect of the post-TOV terms that enter 
in the pressure equation. Right panels: Effect of the post-TOV terms that enter in the mass equation. The circles represent 
contours of fixed relative deviation from GR. 


to the mass-radius curves. 

As a rule of thumb, the corrections to the pressure 
equation lead to markedly weaker mass-radius modifica¬ 
tions than the corrections to the mass equation, for the 
same magnitude of and /i^. The effective-metric for¬ 
mulation of the post-TOV formalism suggests a simple 
qualitative explanation of this observation. The mass 


corrections M 2 change both the effective EOS and the 
strength of gravity, as measured by v{r)^ while the pres¬ 
sure corrections V 2 are only associated with a change in 
the strength of gravity [cf. Eqs. (73) and (76)], and it 
is well known that changes in the EOS outweigh gravity 
modifications in terms of their effect on the mass-radius 
relation. 
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(/^l ^ 0 

(7r2,7r3,7r4) > 0 

(A‘3.A‘4)< 0 

(7r2.7r3,7r4) < 0 

(7ri,A‘3.A‘4>A‘3)> 0 

(a‘i,A‘2)< 0 


{R—Rq^)/Rq^ 


FIG. 9. Directions of the post-TOV induced deviations. This 
schematic diagram shows which sign of individual post-TOV 
parameters produces smaller or larger masses/radii with re¬ 
spect to GR, cf. Fig. 8. 


A notable exception is the single-member family FI, 
for which the pressure correction term dominates over 
its mass counterpart. In fact, the FI pressure term leads 
to the largest mass-radius changes, as evidenced by the 
TTi values used in Fig. 2. It is not too difficult to ex¬ 
plain why this happens: near the stellar surface, where all 
three fluid parameters p, p, 11 are close to zero, the FI cor¬ 
rection terms remain finite and dominate over all other 
terms in the post-TOV equations (this can be clearly seen 
in Fig. 5), thus taking control of the pressure and mass 
derivatives. 

Another noteworthy point is that, when considering 
individual post-TOV terms, it is not always possible to 
integrate the equations for both positive and negative 
values of the corresponding coefficient. This is the case 
for family F5 in Fig. 2, where the integration fails for 
ps < 0. We have found that this is caused by an unphys¬ 
ical negative slope dm/dr near the origin. 

The remarkable self-similarity in the radial profiles of 
2PN terms belonging to the same family has been il¬ 
lustrated in Fig. 3 (see Section IIIB). With hindsight, 
this property should not come as a total surprise, given 
the approximate correlations among the fluid variables: 
m ~ pr^, n ~ p/p. 

The emergence of the same self-similarity in the mass- 
radius curves is something far less anticipated and even 
more striking. This property, which has allowed us to 
formulate a practical and versatile set of post-TOV equa¬ 
tions, is illustrated in Figs. 6 and 7, where we show mass- 
radius results for each 2PN family, considering both the 
pressure and the mass equation and for the same APR 
EOS stellar model as in Fig. 2. Each panel is devoted to 
a particular family, and it shows the mass-radius curves 
resulting from the integration of the post-TOV equations 


when various terms from Table I are included as correc¬ 
tions (notice that FI is missing from these plots for the 
obvious reason that it consists of only one post-TOV cor¬ 
rection). 

In all cases considered, the terms of the same family 
are found to cause nearly identical mass-radius changes 
by a suitable rescaling of the relevant coefficient or pi. 
This behavior is most striking for family F4, where dif¬ 
ferent post-TOV corrections in the mass equations lead 
to the same characteristic back-bending behavior in the 
mass-radius curve. The only notable exception to this re¬ 
markable scaling property is the (0, 0, 1) member of the 
F3 family, proportional to Ilm/r, which can be rescaled 
to agree with other members of the family at high den¬ 
sities but partially fails to capture the behavior of the 
mass-radius curve at low densities. This partial symme¬ 
try breaking can be understood by looking at the leftmost 
panel in the second row of Fig. 3: the behavior of this 
term near the surface is not as smooth as for other mem¬ 
bers of this family. In our opinion this does not warrant 
extensions of the formalism to include another family, 
but this is definitely a possibility that could be consid¬ 
ered in the future, given the approximate nature of the 
self-similarity argument. 

Another important aspect of the post-TOV results is 
their “directionality” in the mass-radius plane, in the 
sense that a given correction term could affect the mass 
more than the radius, or vice versa. This kind of in¬ 
formation cannot be easily extracted from a traditional 
mass-radius plot such as Fig. 2, but becomes very visible 
if we display the same results in terms of the fractional 
changes SM/Mq^ = {M - Mgr)/Mgr and SR/Rqr = 
{R — Rgr)/Rgr from the corresponding GR values. 

“Dart-board” plots of these fractional changes are 
shown in Fig. 8. The aforementioned directionality of the 
various post-TOV corrections is clearly visible in this fig¬ 
ure. Individual correction terms are seen to drive nearly 
linear departures (at least up to a ~ 10% level) from the 
center of the “board.” Moreover, certain terms are mu¬ 
tually (nearly) orthogonal, although not aligned with the 
mass or radius axis. In some cases this happens between 
the pressure and mass terms of the same family, e.g. fam¬ 
ily F2. In general, the departures from the GR model are 
more isotropically scattered when caused by the correc¬ 
tions AI2 in the dm/dr equation, whereas the pressure 
corrections V 2 are clearly more concentrated near the di¬ 
rection of the mass axis. This behavior fits nicely with 
the effective-EOS interpretation of how AI2 and V 2 cor¬ 
rections change the mass-radius diagram. As expected, 
M .2 corrections affect the stiffness of the effective EOS 
with significant effects on the radius, while V 2 corrections 
change the strength of gravity, and this mostly affects 
how much mass a particular model can support. 

These trends remain unchanged as the central energy 
density (and the stellar mass) increases (see bottom pan¬ 
els of Fig. 8). The pressure correction term associated 
with TTi (family FI) provides the exception to the rule: 
a sequence of tti > 0 values leads to a non-linear tra- 
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FIG. 10. Deviations induced by the post-TOV parameters on the stellar compactness. Here we consider the influence of the 
post-TOV parameters on the compactness C = M/R of neutron stars. Deviations from GR are calculated assuming the same 
APR EOS models as in the top and bottom rows of Fig. 8. Left panels: Effect of the post-TOV terms appearing in the pressure 
equation. Right panels: Effect of the post-TOV terms appearing in the mass equation. 


jectory, with initially just the radius decreasing and then 
followed by a comparable fractional decrease in the mass. 
Negative values of tti are not shown because they lead to 
unphysical models where in the outer low-density layers 
of the star dp/dr becomes nearly zero but never negative, 
thus preventing us from finding the exact location of the 
surface (as we have pointed out earlier in this section, 
this behavior is related to the non-zero value of the FI 
term at the surface). 


Fig. 9 provides a schematic chart of the correlation be¬ 
tween the sign of the tt^, pi coefficients and the sign of 
the associated variations 5M, SR. Interestingly, the tt^ 
terms are limited to just two of the four possible quad¬ 
rants (note the anti-correlation between the signs of tti 
and the other tt^). This translates to variations that si¬ 
multaneously make the star bigger (smaller) and heavier 
(lighter), i.e. SR > 0, SM > 0 (or SR < 0, SM < 0). 
In contrast, the pi terms occupy all four quadrants, with 
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the F3, F4, F5 families leading to 5R5M > 0 varia¬ 
tions, and FI, F2 giving rise to the opposite arrangement, 
5R5M < 0. 

The linear patterns of Fig. 8 suggest that the mass and 
radius variations, for a given = {tTj ^ tti, /Xi}, obey the 
empirical relations. 


5M 

Mqb. 




6R 

Rgr 




(77) 


where the structure parameters Km^ Kr are functions of 
the EOS and of Cc, but they are independent of ai. Given 
the nonlinear character of the post-TOV equations, this 
conclusion is clearly nontrivial. We can recast this result 
in terms of the variation of the stellar compactness C = 
MjR, 


— « a, {Km - Kr) . (78) 

t-GR 

This almost linear 5C{ai) dependence^ can indeed be seen 
in the numerical results shown in Fig. 10, where we con¬ 
sider the same stellar models as in Fig. 8. 

The results presented in this section provide a wealth of 
information on the character of the post-TOV corrections 
on stellar structure. It is likely that a more systematic 
study of the self-similar F-families will reveal additional 
layers of information and provide clues as to why the 2PN 
terms change the bulk properties of the star the way they 
do, as a function of the central density. Such a study is 
beyond the scope of this paper but provides an attractive 
subject for future work. 


V. CONCLUSIONS AND OUTLOOK 

This paper is a first step towards establishing a 
parametrized perturbative framework that should, at 
least in principle, encompass all modifications to the bulk 
properties of neutron stars induced by modified theo¬ 
ries of gravity. As in the original formulation of the 
PPN formalism, along the way we were forced to make 
some reasonable simplifying assumptions in order to re¬ 
duce the complexity (and increase the practicality) of our 
parametrization. These reasonable assumptions may well 
fail to match the well-known creativity of theorists, and 


^ It is interesting to note that the qualitative effect of the post- 
TOV terms in the pressure equation can be understood by 
analogy with the case of anisotropic stars in GR. The post- 
TOV pressure equation takes the form dp/dr = {dp/dr)Q^ — 
prmzifi{r)/r'^, with fi{r) > 0, whereas anisotropic stars obey 
dpr/dr = {dpr/dr)Q^ — ‘lair, with a = pr — Pq being the dif¬ 
ference between the radial and tangential pressure. These two 
expressions can be matched if 2ra = prmTifi{r). The compact¬ 
ness of anisotropic stars is known to decrease (increase) when a 
increases (decreases) [61]. This conclusion is in good qualitative 
agreement with the results shown in the left panel of Fig. iO. 


it will be interesting to see how the formalism can be 
extended and improved. 

In a follow-up paper we will use our basic post-TOV 
equations to recover stellar structure calculations in some 
popular theories of gravity, such as those shown in 
Fig. I. It is particularly interesting to compare the for¬ 
malism against theories that violate some of our basic 
assumptions, such as scalar-tensor gravity with sponta¬ 
neous scalarization (which introduces intrinsically non- 
perturbative effects [32]) or Eddington-inspired-Born- 
Infeld gravity, with its lack of a Newtonian limit and 
its unorthodox dependence on the stress-energy ten¬ 
sor [27, 62]. 

We have already obtained some interesting results in 
this context: for example, our conclusion that the 2PN 
post-TOV equations are equivalent to an effective mod¬ 
ified perfect-fluid EOS (see Section IIIC) has an inter¬ 
esting parallel with the results by Delsate et al. [62], 
who reached a similar conclusion for Eddington-inspired- 
Born-Infeld gravity. We are currently extending the “ef¬ 
fective metric” formalism developed in this paper to the 
exterior spacetime of compact stars [60]. This is neces¬ 
sary to compute physical observables such as the gravi¬ 
tational redshift of surface atomic lines, the touchdown 
luminosity of a radius expansion burst and the appar¬ 
ent surface area of neutron stars [63], and it is possible 
that the combination of multiple observables may lift the 
EOS/gravity degeneracy. 

There are several interesting extensions of our work 
that should be addressed in the future. The most obvi¬ 
ous one is to assess whether post-TOV parameters can 
indeed reproduce the mass-radius curve in various classes 
of alternative theories, and whether the post-TOV pa¬ 
rameters encode specific information on the physical pa¬ 
rameters underlying specific theories. This study will 
hopefully lead to a better understanding of the general¬ 
ity of the EOS/gravity-theory degeneracy. 

From a data analysis point of view, it is important 
to understand whether physical measurements of masses 
and radii (or perhaps more realistically, measurements of 
masses and surface redshifts/stellar compactnesses) can 
lead to constraints on the post-TOV parameters under 
specific assumptions on the high-density EOS. The an¬ 
swer to this question obviously depends on the relative 
magnitude of modified gravity effects and EOS uncertain¬ 
ties. It will be interesting to quantify what uncertainties 
in the EOS are acceptable if we want to experimentally 
constrain post-TOV parameters at meaningful levels. 

Other obvious extensions are (i) the generalization of 
the post-TOV framework to slowly and possibly fast ro¬ 
tating relativistic stars, and (ii) stability investigations 
within the post-TOV framework. We hope that our 
work will stimulate further activity in this field. Sta¬ 
bility studies in a post-TOV context may reveal that 
certain generic features of modified gravity lead to in¬ 
stabilities even for nonrotating stars, possibly excluding 
whole classes of modified gravity theories. 

Last but not least, we would like to point out that our 
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post-TOV toolkit is not (nor was it designed to be) a self- 
consistent PN expansion, but rather a phenomenological 
parametrization of the leading-order (unconstrained) de¬ 
viations from GR. A systematic and self-consistent PPN 
expansion extending the PN stellar structure works cited 
in the introduction [40-43] is an interesting but quite dis¬ 
tinct area of investigation that should also be pursued in 
the future. 
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Appendix A: Dimensional analysis of 
post-Newtonian terms 

In this appendix we develop an algorithm for con¬ 
structing PN terms using dimensional analysis tech¬ 
niques. Barring PN terms involving the three poten¬ 
tials U, E, U, the available parameters for generating 
PN terms are {p, p, to, r, H}. From these quantities plus 
the gravitational constant G and the speed of light c we 
can build the dimensionless combination®: 

A = (Al) 

for a suitable choice of integers a,/3,7, i5, k, A (these are 
not to be confused with the PPN parameters of Sec¬ 
tion II). Since H is already dimensionless, there is no a 
priori dimensional restriction on 9 (apart from one com¬ 
ing from the PN order of A) and therefore that factor 
can be omitted in the dimensional analysis. Using the 
scalings 

TO ^ pr®, (A2) 

r 


® Note that this combination is oblivious to the presence of dimen¬ 
sional coupling constants that might appear in modified theories 
of gravity. 


we obtain the following form for A in terms of mass, 
length and time dimensions: 

A~ ^]\^^ja+l3+'r-K^j^j-a+S-3/3+\+3K^rpj-2a-\-2K (A3) 

Since A is required to be dimensionless, we have the three 
algebraic relations: 

A = —2(a-l-K), K = a + /3 + j, (A4) 

and 

-a + S-3f3 + X + 3K = 0. (A5) 

The first two relations simply express A and k in terms 
of the other parameters. Using them in (A5) we obtain 

j + S = 2{a + (i), (A6) 

which represents the true dimensional degree of freedom. 
It is straightforward (if tedious) to verify that all PN 
terms appearing in the PPN equations of Section II are 
consistent with (A6). 

All a < 0 terms are divergent at the surface and need 
not be considered. As we shall shortly see, all terms 
with a > 4 are divergent at r = 0 in both structure 
equations, and therefore should be discarded. The a = 
3 terms are singular in the dp/dr equation and can be 
discarded by the same argument; a = 3 terms are regular 
in the dm/dr equation, but they are always dominated 
in magnitude by the a < 3 terms, and therefore will 
not be presented in detail here. Therefore our strategy 
hereafter is to focus on the particular cases 0 = 0 (no 
pressure dependence) and a = 1, 2 (linear and quadratic 
scaling with the pressure). 

1. Terms with a = 0. 

Starting with the a = 0 case we have 

7 -b <5 = 2/3. (A7) 

The resulting form of A in geometric units is 

(AS) 

Formally, this combination is of order (jn/r)^'^"^. There¬ 
fore, we can generate A^-PN terms if /3 -I- 7 = A. These 
are of the form 

(A9) 

with /3 = 0, ±I,±2,.... For instance, the first few IPN 
and 2PN terms of this series are (we start from /3 = — I 
for reasons explained below): 


Ai(-l)^ 


Ai(0)^ 

Ai(l)~r2p, (AIO) 

r 


3 


2 

A2(-1) - 

m 

e: 5 

r^p 

A2(0) ^ 

771 

" ^ 2 ( 1 ) ~ rpm. (All) 
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2. Terms with a = 1. 

The a = 1 group of terms can be obtained with the 
same procedure. We have 

7 + 5 = 2(1 + /3), (A12) 

and this leads to terms of the form 

. (A13) 

Since r^p is a 2PN term, the resulting A^-PN combination 
should take the form: 

An{13) ~ . (A14) 

The first few IPN and 2PN terms generated from this 
expression are: 

Ai(-l)^^, Ai(0)~^, (A15) 

A2(-1) ~ , A2(0) ~ rV A2(1) ~ (A16) 

pr m 

3. Terms with a = 2. 


Finally, we consider the a = 2 terms. The correspond¬ 
ing Ajv combination is. 


^Ni(5) = {r'^p)^ir^p)^ 

(")"■ 

- 4-/3 

5 

(AI7) 

and from this we have: 




4 2 

r V 

Ai(-I)~^, 

m^p 

Ai(0) ^ 

j.7p2 

m^ ’ 

(AI8) 

j,3p2 

A2(-I)^^, 

pm 

A2(0) - 

rSp2 

m^ 

(AI9) 


4. Generic A-PN order terms and constraints. 

It is now not too difficult to see that a A-PN order term 
with an arbitrary p°‘ scaling and with 11 re-introduced is 
given by the universal formula, 

AAr(a,/3,6») - n®(r2p)“(rV)^ j . (A20) 

As discussed in Section IIIB, different threads of reason¬ 
ing lead to the constraint /3 > — 1 The first one has to 
do with avoiding a divergence at the stellar surface (this 
already has allowed us to filter out all a < 0 terms). An 
inspection of the two stellar structure equations reveals 
that terms with a = 6 = Q should scale as 

pAjv(0,/3,0)^pi+^, (A21) 


in the vicinity of the surface, and therefore we ought to 
take /3 > — 1 in order to avoid a surface singularity. This 
argument still allows for /3 < — 1 values in the A^ terms 
with a, 6* > 0, since these terms have a smoother profile 
as a result of the vanishing of p and 11 at the surface. 

The second thread is no more than a heuristic ar¬ 
gument and has to do with the expectation that for a 
broad family of gravity theories the solution for the met¬ 
ric (and its derivatives) should scale as ~ (e -I- rp)” = 
/9”(1 -I- n -f rp/p)" with the fluid parameters (where r 
and n are 0{1) numbers). From this it follows that neg¬ 
ative powers of p will come in the form of dimensionless 
PN terms ^ p'^~^{p/p)^, where k = n, n — 1,... (note 
that a factor p has been absorbed by the Newtonian pref¬ 
actor in the structure equations). As a consequence, 
is the only possible negative power in a PN expansion. 
Obviously, this argument automatically takes care of the 
regularity of any (5,9) term at the surface. 

The exclusion of all a > 4 terms comes about as a 
consequence of regularity at the stellar center. Near 
the origin (where p, p, 11 take finite non-zero values) a 
Aff{a,(3,6) term behaves as 

An^t ^ 0) ~ r^(N-a-e)^ (^22) 


pm 

^2{N-a-e) + l 

(A23) 

' r^pA/v - 

^ ^2(N-a-e+l)^ 

(A24) 


The corresponding terms in the stellar structure equa¬ 
tions will behave as 


dr 

dm 
dr 

and therefore regularity at the center dictates the follow¬ 
ing limits for each equation: 

y : 0<a<N-e, (A25) 

dr 

fjrn 

V: O<a<A+l-0. (A26) 

dr 

We can also see that these conditions entail the following 
limits for 9: 


0<e<«, 

dr 


dm 

— : Q<9<N+1. (A27) 

dr 


0 <6»<2, 0<a<2-6», (A28) 


For the particular case of 2PN order terms we then have: 
dp 
dr 
dm 
dr 

which shows that all a > 4 terms are to be excluded and 
that a = 3 terms can only appear in the mass equation. 


0<6»<3, 0<a<3-i 


(A29) 


Appendix B: The Newtonian and relativistic 
Lane-Emden equations 

In this appendix we review the nonrelativistic and rel¬ 
ativistic Lane-Emden equations. The former equation is 
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classic textbook material (see e.g. [64]) and therefore is 
just sketched here. The somewhat less familiar relativis¬ 
tic extension was developed by Tooper [65, 66] and is 
discussed in a bit more detail below. Our definition for 
the polytropic EOS, i.e. p = is the same as 

the one adopted in [66] but is different to the one used in 
Tooper’s earlier paper [65], i.e. p = This sub¬ 

tle difference, combined with the choice between pdPc 
or Pcjic (the “c” index refers to the stellar center) for 
the scale of the system, leads to slightly different Lane- 
Emden equations. 


1. The Newtonian Lane-Emden equation 


In Newtonian gravity, one can express the hydrostatic 
equilibrium equation for spherical non-rotating stars in 
terms of dimensionless parameters for the pressure, the 
density and the radial coordinate. If the EOS is poly¬ 
tropic (i.e., according to our definition, p = 
the equations governing the dimensionless quantities are 
scale-invariant, depending only on the polytropic index 
n. By writing the density and the pressure as 


0" = ^, p = Kpl+^^^e^+\ (BI) 

Pc 

and introducing the dimensionless radial coordinate 


r = aC, 


a = 


(n -|- 1)K 
47rG 


1/2 




(B2) 


the Newtonian stellar structure equations 


lead to 


dp 

dr 

dmN 

dr 


Gmjq 

J.2 

47rr^p, 


P, 


dd^ V di) 


0 " 


(B3) 

(B4) 


(B5) 


This is the famous Lane-Emden equation, and its scale- 
invariant solutions describe all possible fluid configura¬ 
tions in terms of the single parameter n. 


2. The relativistic Lane-Emden equations 

Generalizing the Lane-Emden formalism to GR is a 
straightforward task, but this comes at the price of los¬ 
ing the scale-invariance property of the Newtonian treat¬ 
ment. In relativity we can dehne the polytropic EOS in 
the same way as before, where p is the baryonic rest mass 
density. The polytropic exponent is defined as 


Then the energy density e and the internal energy 11 are 
given by 


e = p + np, (B7) 

which implies 

U = n-. (B8) 

P 

This observation was used in the argument leading to 
Eq. (53). 

We can now introduce the relativistic version of the 
Lane-Emden equations. In analogy with the Newtonian 
case we define p = Pc0", r = and p = 

The ratio between the central pressure and the central 
energy density 



H-l/n 

Kpc 

I IV 1 + 1/' 

Pc + nKpc 


(B9) 


is a convenient measure of the importance of relativistic 
effects in the system. Note that our definition deviates 
from Tooper’s [66], who prefers to use the ratio pd Pc- 
The energy density is then 

e = PcO'^ + nKp]+^/'^e'^+^ 

= ec[I-fnA(0-I)]0". (BIO) 

We now want to derive a dimensionless form of the 
TOV equations (7). The definition of the mass function 
mx implies 


= 47reca^[I -I- nX{9 — I)]0"^^. 

a? 

In terms of the dimensionless mass 

mx 


this becomes 


dm 


47r[l-hnA(0- 1)]0”^2. 


(BIl) 


(BI2) 


(BI3) 


From the TOV equation for the pressure we similarly 
obtain, after some manipulations. 


d9 m 


l + (n + l) — 


nX 


X I 1 + A 


4^^35,n+l 


m 


TD 

l-2{n + l)Xj 


(BI4) 


In the present case the characteristic length scale is 

1/2 


a = 


{n + l)Kpd^^^^d-nXy 


(BI5) 


1 p dp e+ p dp 
n p dp p de 


(B6) 


At this point we would like to define dimensionless 
quantities that come from the relativistic Lane-Emden 
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equations. The central baryonic rest-mass density is re¬ 
lated to A as [see Eq. (B9)]: 

p, = K-^r, e=-^. (B16) 

The factor K~^ has units of mass density (or inverse 
square length in geometrical units), therefore the dimen¬ 
sionless rest-mass density is 


p = pK'^ = re^. (B17) 

Similarly, the length scale a takes the form 

a = + 1)^^-” (1 - nX) , (B18) 

where has dimensions of length. The dimensionless 
radius is defined as 

f = (1 - nA) (B19) 

The remaining dimensionless parameters are 

e = eif” = [l + nA(0-l)]r, (B20) 


p = 

= (1 — nA) 


1 — nA 


TO, 


(B21) 

(B22) 


p = = r+i6»”+\ (B23) 

n = n^=n£6 ». (B24) 

P 

All of the above dimensionless profiles are functions of 
n and A. At variance with the Newtonian treatment, the 
relativistic Lane-Emden formalism does not allow for a 
simple algebraic mass-radius relation M{R). This is also 
related to the fact that the system is not scale-invariant, 
due to the presence of A in the equations. 


Appendix C: The PPN Potentials 



FIG. 11. Integral PN potentials and IPN terms. The radial 
profiles of the integral potentials U, E/m-[q and are well 

fitted by linear functions of the non-integral potentials p/p 
and r^p/mjq. In all plots, the radial coordinate is normalized 
to the stellar radius R. 


where all right-hand side quantities are computed in 
Newtonian theory. In Eqs. (Cl), TON('r) denotes the New¬ 
tonian mass function 


ton(?') = An dr' pr'"^ 

Jo 


r 

dTTTOb / dr'n\yr'^^ 

Jo 


(C2) 


The goal of this appendix is to study the behavior of 
the potentials U, E and 17 appearing in the PPN stel¬ 
lar structure equations (9), first derived by Ciufolini and 
Ruffini [47]. By means of a mass function redefinition (see 
Section II) these potentials can be eliminated at IPN or¬ 
der, but they could still appear at 2PN order and higher. 

Given the 2PN precision of our calculations we can 
write these potentials as: 

U{r) = - r dr'^^ + UiO), (Cla) 

Jo ^ 

E{r) = An I dr'r'^pH, (Clb) 

Jo 

^{r) =—An f dr'r'pm^, (Clc) 

Jo 


where rib is the baryon number density. The integral 
quantities U, E and 17 represent the system’s gravita¬ 
tional potential energy, internal energy, and gravitational 
potential energy respectively [34] . They appear as dimen¬ 
sionless PN terms in the form of reduced potentials: U, 
E/mys, f7/TON [see Eqs. (9)]. 

The radial profiles of the three potentials inside the 
star can be determined by first integrating the Newtonian 
hydrostatic equilibrium equations, Eqs. (B3) and (B4), to 
find Ton and p as functions of r. Using realistic EOS data 
tables for p{p) we can subsequently compute the internal 
density per unit mass 11 (p) and the mass density p{p) = 
Wbiib(p), and then numerically evaluate the potentials 
inside the star by integration. 

Some insight into the nature of these potentials can be 
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obtained by rewriting Eqs. (Cl) in the form 


u 

= — +Att [ 

dr'r' p. 

(C3a) 


r 

J r 



E 

1 

pr 

dn 



= n- 

1 dr' 

(C3b) 
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_ Wn 
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Ton 

2r 

2ton 

Jo \ r' J 



= 4.^- 

127r 

[ dr' p. 

(C3c) 


TON 

Ton 

Jo 



Note that the integration constant for U has been fixed 
by requiring U{R) = M/R at the stellar surface, while 
those for E and il have been set to zero in order to have 
regularity of if /ton and JI/ton at r = 0. The values of 
the potentials at the stellar center are: 

t/(0)=47r / drrp, -(0) = 0, -(0) = lie. 

Jo Ton 

(C4) 

From Eqs. (C3) we can see that E/m-^ and H/ton are 
(partially) expressed in terms of the non-integral IPN 
terms 


WN „ ^ 
r ’ ’ mN 


(C5) 


This suggests the possibility that the behavior of all 
three potentials could be captured by linear combina¬ 
tions of non-integral IPN terms. If true, this would mean 
that any 2PN term involving [/, E/ton or H/ton is effec¬ 
tively accounted for by the presence of the other terms 
in the post-TOV formulas. For instance, this idea can be 
demonstrated for U and for the special case of a poly¬ 
tropic system. Starting from (Cla) and expressing ton 


in terms of dp/dr, after an integration by parts and use 
of (B6) we arrive at 


U={n + l)(--—]+U{0). (C6) 

\P PeJ 

We know that for a polytrope 11 = np/p, which means 
that we can also write 

(^ + 1) (n-Ee)-H/(0). (C7) 

n 

For a polytropic model, therefore, U can be written ex¬ 
actly as a linear function of p/p or 11. 

We have verified that U, E/m-^ and H/ton can be 
approximated by similar linear functions for the case of 
realistic EOSs. As an illustration, in Fig. 11 we con¬ 
sider a stellar model built using the APR EOS with cen¬ 
tral mass density of 0.58 x 10^® g/cm^, Newtonian mass 
Ton = I. 5OM0 and radius R — 14.8 km. For this model 
we plot the radial prohles of U (top panel), E/ton (mid¬ 
dle panel) and 0 /ton (bottom panel). The figure shows 
that the profiles of the three potentials can be accurately 
reproduced by linear combinations of the IPN terms in 
Eq. (C5), and that U is reasonably well fit by a linear 
function of p/p, as suggested by (C6). This latter fit 
breaks down near the surface, but with a different com¬ 
bination of IPN terms (namely, 11 and r^p/m^) one can 
produce a near-perfect fit. 

In conclusion, the addition of the integral potentials 
U, E/rriN and H/ton in the 2PN terms is unnecessary 
because their behavior can be captured by linear combi¬ 
nations of the non-integral PN terms which are already 
included in the post-TOV equations (1). 
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